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ABSTRACT 

We present results from an analysis of Far Ultraviolet Spectroscopic Explorer (FUSE) spectroscopy of 
the z em — 0.57 quasar PKS0405-123. We focus on the intervening metal-line systems identified along 
the sightline and investigate their ionization mechanism, ionization state, and chemical abundances. In- 
cluding HST/STIS spectroscopy, we survey the entire sightline and identify six OVI absorbers to a 3<r 
equivalent width (EW) limit of 60mA. This implies an incidence dN/dz — 16^g consistent with previous 
O VI studies. In half of the O VI systems we report positive detections of C III suggesting the gas is pre- 
dominantly photoionized, has multiple ionization phases, or is in a non-equilibrium state. This contrasts 
with the general description of the warm-hot intergalactic medium (WHIM) as described by numerical 
simulations where the gas is predominantly in collisional ionization equilibrium. An appreciable fraction 
of OVI absorbers may therefore have a different origin. We have also searched the sightline for the 
NeVIII doublet over the redshift range 0.2 < z < z em which is a better probe of the WHIM gas at 
T > 10 6 K. We find no positive detections to an EW limit of 80mA giving dN/dz < 40 at 95% c.l. 

The photoionized metal-line systems exhibit a correlation between the ionization parameter (U = 
&/criH with <!> the flux of hydrogen ionizing photons) and HI column density for iV(HI) = 10 14 — 
10 16 cm -2 . Both the slope and normalization of this correlation match the prediction inferred from the 
results of Dave and Tripp for the low z Lya forest. In turn, the data show a tentative, unexpected 
result: five out of the six photoionized metal-line systems show a total hydrogen column density within 
a factor of 2 of 10 'cm -2 . Finally, the median metallicity [M/H] of twelve z ~ 0.3 absorbers with 
iV(HI) > 10 14 cm~ 2 is [M/H] > —1.5 with large scatter. This significantly exceeds the median metallicity 
of CIV and OVI systems at z ~ 3. and therefore requires enrichment of the intergalactic medium over 
the past ss 10 Gyr. 



Subject headings: Intergalactic medium 
1. INTRODUCTION 

Quasar absorption line studies examine the gas within 
galaxies and the intergalactic medium (IGM) between them. 
These studies include analyze of the interstellar medium 
of young galaxies at z > 2 (e.g. the damped Lya systems) 
as well as the lines that comprise the Lya forest. One of 
the most notable successes of cold dark matter cosmolo- 
gies at z ~ 3 is the general agreement between observa- 
tions of the Lya forest and cosmological simulations both 
in terms of the column density and Doppler width dis- 
tributions (e.g. Miralda-Escude et al. 1996; Zhang et al. 
1997). Meanwhile, analyses of metal- line systems describe 
the enrichment history of the universe (Pettini et al. 1997; 
Prochaska et al. 2003), test processes of nucleosynthesis 
(e.g. Lu et al. 1996; Prochaska, Howk, & Wolfe 2003), 
and examine the role and characteristics of winds (Aguirre 
et al. 2001) and Population III enrichment (Gnedin & Os- 
triker 1997; Wasserburg & Qian 2000). 

The 'fluctuating Gunn-Peterson' paradigm inferred from 
numerical simulations describes the majority of the ob- 
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served IGM as large-scale overdensities with little corre- 
spondence to individual galaxies (e.g Gnedin & Hui 1998). 
Although this model reasonably reproduces the number 
density and gross kinematics of the IGM, it has not been 
extensively tested by observations. This is especially true 
at low redshift where there are very few sightlines with 
high resolution, high quality spectra. For this reason, al- 
ternative models (e.g. absorption by gas in galactic halos 
or low surface brightness galaxies; Linder 2000; Chen, 
Lanzetta, & Webb 2001; Manning 2003) are generally 
as successful at explaining the low z IGM as the CDM 
paradigm (Dave & Tripp 2001). 

We are pursuing an observing program aimed at ad- 
dressing key questions related to the low z IGM. We have 
chosen to study the IGM along the sightlines to z > 0.1 
AGN with relatively high signal-to-noise (S/N) far-UV 
FUSE spectroscopy and, in many cases, HST/STIS echcllc 
observations. This UV spectroscopy is complemented by 
large field of view (> 20' diameter) galaxy surveys com- 
piled with the WFCCD spectrograph at Las Campanas 
Observatory. The principal focus of our program is to iden- 
tify the location and physical characteristics of low redshift 
OVI absorbers, whose hot gas content is believed to be 
a significant reservoir of baryons (Cen & Ostriker 1999; 
Tripp et al. 2000; Dave et al. 2001; Fang & Bryan 2001). 
The properties in which we are interested include the ion- 
ization state, metallicity, and absorber/galaxy cross-correlation 
function for OVI gas. Studies of the OVI absorbers cur- 
rently provide the most efficient means of studying the 
relatively hot T k, 10 5 ~ 7 K, low density gas termed the 
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warm- hot intergalactic medium (WHIM). The WHIM is 
predicted to contain a large fraction of the baryons at low 
z (Dave et al. 2001), but its existence is difficult to ver- 
ify due to its low density. In addition to studies on the 
origin of the WHIM, the combined datasets of high qual- 
ity UV spectroscopy and deep, large field-of-view galaxy 
surveys will examine several aspects of the IGM/galaxy 
connection. These issues are intimately related to our 
understanding of gas and galaxies in the local universe. 
Through comparisons with hydrodynamic numerical sim- 
ulations, we will test and constrain the CDM paradigm 
for the IGM and examine processes of galaxy feedback 
and chemical enrichment. 

This paper focuses on the metal-line absorbers identi- 
fied in the ultraviolet spectroscopy of PKS0405-123. We 
have focused first on the PKS0405-123 sightline for the 
following reasons: (i) It exhibits one of the brightest UV 
fluxes for a z > 0.5 quasar. This allowed the STIS sci- 
ence team to obtain high quality UV spectroscopy from 
A « 1200 — 1700 A and our team to acquire modest S/N 
FUV spectroscopy from A « 900 - 11 70 A with the FUSE 
Observatory, (ii) It exhibits a partial Lyman limit system 
(i.e. an absorber with r < 1 at \ res t = 912A) at z = 0.167 
(Chen & Prochaska 2000). In this paper we report on our 
analysis of the UV spectroscopy and examine the metal- 
line systems identified along this sightline. and (iii) We 
have composed a « 1D° galaxy survey reaching R < 20 
with w 500 galaxies at z < z em . Concurrent and forthcom- 
ing papers on this sightline will (a) report on the physical 
nature of the absorber identified at z = 0.4951 (Howk et 
al. 2004), (b) describe the galaxy survey and the connec- 
tion between galaxies and the metal-line systems discussed 
here (Prochaska et al. 2004); and (c) present an analysis 
of the Lya/galaxy cross-correlation function (Chen et al. 
2004) . Future work will analyze the IGM and the galaxies 
associated with it for « 10 low redshift fields. 

The paper is organized as follows: in § 2 we present the 
FUSE observations and a line list of significant absorbers 
identified along the sightline; in § 3 we describe our ap- 
proach to ionization modeling; we analyze the ionization 
state, physical conditions and elemental abundances of the 
metal-line systems in § 4; § 5 lists the absorbers along the 
sightline with iV(HI) > 10 14 cm~ 2 which do not display 
metal transitions; we discuss the implications of our re- 
sults for the WHIM, chemical enrichment, and the general 
IGM in § 6; and § 7 gives a brief summary. 



program was principally motivated by the identification of 
a strong metal-line system at z = 0.167 in the HST+STIS /E140M 
data of PKS0405-123 obtained by the STIS science team 
(Williger et al. 2004). In Chen & Prochaska (2000), we 
analyzed this absorption system and argued that it is a 
partial Lyman limit system with roughly solar metallicity. 
We then obtained far-UV spectroscopy with the FUSE ob- 
servatory to further constrain the physical characteristics 
of this absorber. Table 1 summarizes the UV spectroscopy 
acquired with FUSE and HST/STIS of PKS0405-123 and 
estimates the signal-to- noise (S/N) per resolution element. 



Table 2 
ATOMIC DATA 



Transition 


A 


/ 


Rcf 


HI 972 


972.5368 


0.029 


1 


OI 988a 


988.5778 


0.0005146 


1 


OI 988b 


988.6549 


0.007712 


1 


OI 988 


988.7734 


0.04318 


1 


NIII 989 


989.7990 


0.1066 


1 


Sill 989 


989.8731 


0.133 


1 



References. — Key to References — 1: Morton (1991); 2: Howk et al. (2( 
Lawlcr (1999); 7: Fedchak, Wiese, & Lawler (2000); 9: Schectman, Povolny, 
Fcdchak, & Lawler (2001); 13: Bergeson & Lawlcr (1993); 14: Bcrgeson, Mul 
Mullman, & Lawler (1996) 

Note. — [The complete version of this table is in the electronic edition of tl 



Table 3 

ADOPTED SOLAR ABUNDANCES 



Elm 


e(X)- 


Z 


H 


12.00 


1 


C 


8.59 


6 


N 


7.93 


7 


O 


8.74 


8 


Ne 


8.08 


10 


Mg 


7.58 


12 


Al 


6.49 


13 


Si 


7.56 


14 


S 


7.20 


16 


Fe 


7.50 


26 



e(X) is the logarithmic number density of a given element scaled such that e(I 

PKS0405-123 was observed with the FUSE observa- 
tory (Moos et al. 2000; Sahnow et al. 2000) using the 
LWRS apertures with the detectors in TTAG mode dur- 
ing one visit beginning UT 05 October 2001. We used 
ttag_combine to concatenate the photon event lists into a 
single file and then processed this file with the CalFUSE 
a Signal-to-noise per resolution clement for FUSE data and per pixel for HSffiwSise (v2.4). The reduction proceeded without any sig- 
nificant warnings and the standard series of output files 
was created. We binned these data arrays by three 'pix- 
2. OBSERVATIONS, DATA REDUCTION, AND EW ANALYSIS els ' (~ 0.028A per bin) and then traced the continuum of 

PKS0405-123 in each channel separately by fitting Leg- 
PKS0405-123 was observed for 71ks during Cycle 2 of endre polynomials to regions of unabsorbed quasar con- 
thc FUSE mission (Program B087; PI: Prochaska). This tinuum. These normalized spectra form the dataset for 



Table 1 
UV SPECTROSCOPY 



Instr. 


Mode 


Wavelength 


Exp. 


S/N a 


FUSE 


LWRS 


900-1180A 


71ks 


5-15 


HST/STIS 


E140M 


1150-1700A 


27ks 


7 
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Table 4 
EW SUMMARY 



Ion 






^abs 






W 2 


a{W 2 ) 


Wt 


a(W f ) 




(A) 


(A) 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


FUSE 


OIV 787 


931.521 


787.711 


0.18257 


7 


15 


42 


12 


28 


9 


OIV 787 


931.797 


787.711 


0.18292 


124 


17 


131 


11 


128 


11 


MgX 624 


934.389 


624.950 


0.49514 


-2 


13 


-12 


10 


-8 


8 


OV 629 


941.457 


629.730 


0.49502 


183 


12 


145 


12 


164 


9 


SIV 809 


944.893 


809.668 


0.16701 


13 


16 


31 


15 


23 


11 


Oil 834 


973.896 


834.466 


0.16709 


8 


17 


17 


12 


11 


10 


OIII 832 


985.062 


832.927 


0.18265 


34 


15 


5 


11 


15 


9 



Note. — Note that the list is incomplete for wavelengths < 1000A where the data has poor S/N and significant line blending. 

Note. — Columns 4,5 (6,7) refer to the SiClB (SiC2A) channel for 910A< A < 990A, LiFIA (L1F2B) for 995A< A < 1057A, LiFIA for 
1057A< A < 1082A, L1F2A for 1082A< A < HOOA, LiFIB (L1F2A) for 1100A< A < 1180A, and STIS/E140M for A > 1190A. 

Note. — [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 



our analysis. The atomic data considered in this paper is 
tabulated in Table 2 and Table 3 lists our assumed solar 
abundance data compiled by Grevesse & Sauval (1999) 
and revised by Holweger (2001) for C,N, and O. 

Redward of 1000A, we identified all 3tr features in the 
LiFIA and LiF2A channels. There is moderate absorption 
from Galactic molecular hydrogen along this sightline; we 
attribute « 50% of the features at A = 1000 - 1150A 
to H2 or Galactic metal-line transitions. Table 4 lists 
the observed wavelength, transition name, absorption red- 
shift, rest equivalent width W\, and the statistical error 
in W\ from each channel where the transition is observed 
at reasonable S/N. Columns 4,5 (6,7) refer to the SiClB 
(SiC2A) channel for 910A< A < 990A, LiFIA (LiF2B) 
for 995A< A < 1057A, LiFIA for 1057A< A < 1082A, 
LiF2A for 1082A< A < HOOA, LiFIB (L1F2A) for 1100A< 
A < 1180A, and STIS/E140M for A > 1190A (see Chen 
& Prochaska 2000 for details of the STIS dataset and re- 
duction). Note that there are offsets between the wave- 
length solutions for each FUSE channel. We corrected for 
these offsets with a cross-correlation analysis allowing for 
a zero point and a first order term. The latter term, while 
small, is required to match the spectrum over the entire 
spectral range of each channel. This procedure only cor- 
rects for relative errors; a comparison of HI Lyman lines 
and metal lines from the STIS and FUSE spectra indi- 
cates the absolute wavelength calibration is accurate to 
better than 10 km s _1 . The final columns in Table 4 list 
the adopted W\ value and uncertainty corresponding to 
the variance weighted mean for transitions with multiple 
measurements. In general, the values derived from multi- 
ple channels show good agreement. For several transitions 
(e.g. Galactic PII 1152), however, the individual W\ val- 
ues differ by > 3cr statistical significance. In most cases, 
this may be attributed to fixed pattern noise in the detec- 
tors or continuum uncertainties. The transitions showing 
these discrepant cases were either ignored or we adopted 
an uncertainty equal to the spread in the two values. 



X 



o 




Fig. 1. — Predicted ionic column densities for a photoionized 
gas as a function of ionization parameter U . The curves are the 
results of a series of calculations with the Cloudy software pack- 
age (Ferland 2001) assuming (i) a quasar-only extragalactic UV 
background (model Q) at z = 0; (ii) a total HI column density 
Af(HI) = 10 1 ® cm -2 ; and (iii) solar relative abundances scaled to a 
metallicity of 1/10 solar. 



3. IONIZATION MODELING 

3.1. Photoionization 

In this paper, we examine gas likely to be significantly 
ionized. For those systems which show multiple ionic species, 
it is possible to constrain the ionization mechanism (i.e. 
photoionization vs. collisional ionization) and ionization 
state of the gas. In general, we expect the gas to be 
significantly photoionized, either by an extragalactic UV 
background (EUVB) radiation field or a local radiation 
source (e.g. proximity to a starbursting galaxy). For the 
shape of the EUVB field, we consider two models gen- 
erated by Haardt & Madau (2004) using their CUBA 
software package: quasar-only (Q) and quasars+galaxies 
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(QG). These models input the luminosity function of AGN 
and UV bright galaxies for a given redshift and then cal- 
culate the attenuation of the mean intensity by the inter- 
galactic medium. The majority of previous studies have 
adopted model Q for the EUVB field. In our analysis, we 
will consider this model as the default spectrum, and we 
note that ongoing programs studying the UV emission of 
low z galaxies (e.g. GALEX) will provide a more accurate 
assessment of the proper EUVB model. Where relevant, 
we describe the effects of including a softer radiation field 
than model Q. 

In practice we calculate solutions for a plane-parallel 
slab of gas using the Cloudy software package (Ferland 
2001). Because all of the absorbers along the PKS0405- 
123 sightline are optically thin to the radiation field, the 
calculation is simplified and the models are fully described 
by (i) the ionization parameter U = Q/cnn, with $ the 
flux of hydrogen ionizing photons and tih the volume den- 
sity of hydrogen; (ii) the shape of the ionizing radiation 
field; and (iii) the metallicity of the gas [M/H]. The mod- 
els are most sensitive to variations in the ionization pa- 
rameter while the shape of the radiation field and the gas 
metallicity play more modest roles. Throughout the anal- 
ysis we assume constant density (typically an arbitrary 
log Tin = — 1) and constrain the plane-parallel slab to have 
the thickness required to yield the observed iV(HI) value 
given the specific ionization state of the gas. 

We have constructed grids of these photoionization mod- 
els varying U, [M/H], and the EUVB model for a range 
of redshifts and A^(HI) values. An example of one set of 
models is shown in Figure 1 where we have assumed EUVB 
model Q at z = and iV(HI) = 10 16 cm" 2 . The lines trace 
predicted ionic column densities vs. logt/ for a series of 
ions in a gas with 1/10 solar metallicity (i.e. [M/H] = — 1) 
and solar relative abundances (Table 3). Low-ion species 
(e.g. Sill) show less variation at logt/ < —3 because they 
trace the HI gas to first order and the models assume a 
fixed iV(HI) value. In contrast, high-ion species show large 
values only for higher ionization parameters. 

The general effect of adopting a softer radiation field 
(e.g. EUVB model QG) is that one requires a higher ion- 
ization parameter to attain significant column densities of 
CTV,OV, etc. In turn, higher U values imply lower neu- 
tral fractions and therefore larger iV(H) values for a fixed 
(observed) iV(HI) value which - contrary to our initial ex- 
pectation - predicts lower metallicities (~ 0.3 — 0.5 dex 
lower) for absorbers dominated by high-ions. Systems 
characterized by a mixture of low and intermediate ions, 
however, show minimal differences in their elemental abun- 
dances as a function of the EUVB model. 

Our general approach for a given absorber is to (1) con- 
strain U with one or more ionic ratios, ideally of the same 
element (e.g. iV(SiII)/iV(SiIII)); (2) determine if a self- 
consistent solution is allowed with all of the ions included; 
and (3) estimate the ionization parameter, metallicity, rel- 
ative abundances, and physical conditions (e.g. density, 
temperature) of the gas from the best ionization model. 
This technique is dominated by systematic error, e.g. the 
limited constraints on the EUVB model, the simplifying 
assumption of constant density, etc. With current model 
sophistication and typical observational constraints, we do 
not believe one can achieve better than 0.2-0.3 dex preci- 



sion on elemental abundances or 0.2 dex for relative abun- 
dances. Therefore, we take a cautious approach when com- 
paring the observations against the photoionization mod- 
els and adopt conservative errors. 




10 5 10 6 
T(K) 



Fig. 2. — Predicted ionic column densities for a collisionally ion- 
ized gas in equilibrium as a function of gas temperature T. The 
curves were calculated with the Cloudy software package assuming 
(i) no ionizing radiation field; (ii) 7V(HI) = 10 16 cm - 2 ; and (iii) 
solar relative abundances scaled to a metallicity of 1/10 solar. 



3.2. Collisional Ionization 

For gas at relatively high temperatures (T > 5 x 10 4 K), 
collisional ionization can be the dominant ionization pro- 
cess for a number of high-ions observed in quasar absorp- 
tion line systems (e.g. CIV, NV, OVI). In Figure 2 we 
plot abundance curves vs. temperature under the assump- 
tion of collisional ionization equilibrium (CIE; Sutherland 
& Dopita 1993) as calculated by the Cloudy package. 
We have assumed 1/10 solar metallicity and plot the pre- 
dicted column densities for several ions assuming AT(HI) = 
10 16 cm -2 and solar relative abundances (Table 3). Colli- 
sional ionization is important for higher ionization states 
of Si and C at T w 5 x 10 4 - 2 x 10 5 K whereas O and 
N become highly ionized at T w 10 5 — 10 6 K. There are 
several key differences between the curves presented in Fig- 
ures 1 and 2 which direct our discussion on the ionization 
mechanism of the metal-line systems toward PKS0405- 
123. In particular, note that the CIII and OVI ions do 
not coexist with significant values in the collisional ioniza- 
tion solutions yet both ions show large column densities 
for logU ~ —1.5 in the photoionization model. A similar 
difference is noted for NTH. In addition, the relative abun- 
dances of the individual O ions show modest differences 
in the two models, e.g. photoionization tends to predict 
smaller differences between A(OIV) and TV(OVI) than 
collisional ionization. Another key difference between the 
two scenarios is temperature of the gas: CIE models gen- 
erally predict much higher temperature. Therefore, the 
Doppler width (b value) of the absorption can provide a 
discriminant between the two scenarios, especially when 
the neutral hydrogen gas shows b < 30 km s _1 . 
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3.3. Multi-phase Scenarios and Other Ionization States 

If one is limited to the observation of only a handful of 
ions, however, it may be difficult to distinguish between a 
gas with an ionization phase dictated by photoionization 
versus collisional ionization. For example, a CIE model 
with T w 2 x 10 5 K predicts similar relative ionic col- 
umn densities of O VI and N V as a photoionization model 
with logU > — 1. In these cases, degenerate solutions ex- 
ist which limit our constraints on the physical properties 
of the gas. Furthermore, the ionization mechanisms are 
not exclusive. For the absorbers considered in this paper, 
one expects that a radiation field is always present and 
a pure CIE model is somewhat unrealistic (although the 
effects of photoionization may be minor). In the following 
section, our approach will be to assess single-phase ioniza- 
tion models first and only consider multi-phase scenarios 
as required by the observations. 

There are also examples in the literature of ionized gas 
which is not well explained by equilibrium models of col- 
lisional or photionization. The best-known examples of 
non-equilibrium gas involve gas in the interstellar medium 
of galaxies at temperatures where the cooling rate is ex- 
tremely high, including T 10 5 K gas studied through lines 
from the Li-like ions O VI, N V, and C IV. The prevail- 
ing conclusion is that this gas is not in equilibrium and 
arises as a result of conductive heating, radiative shocks, 
and/or turbulent mixing layers (see Spitzer 1996 for a sum- 
mary). Models of these phenomena predict a wide range of 
gas properties and diagnostics dependent on the assumed 
initial conditions and timescales for cooling, heating, and 
mixing. In this paper, we will restrict our analysis equi- 
librium models of collisional and photoionization. In part, 
this is because we have very few measurements of C IV and 
Si IV for these absorbers. More importantly examples of 
conduction heating and radiative shocks are more likely to 
occur in dense, star forming regions; such conditions are 
unlikely for the absorption systems considered here. Fi- 
nally, we will find that the systems are well described by 
a single or two-phase equilibrium model. Nevertheless, it 
is important to stress that our results are sensitive to the 
presumption of equilibrium modeling. 

4. ANALYSIS OF THE METAL-LINE SYSTEMS 

This section presents an analysis of all of the metal-line 
systems identified along the sightline to PKS0405-123. Be- 
cause of the complicated line-spread functions of the spec- 
trographs, we derive the iV(HI) values and Doppler pa- 
rameters for the HI gas using a curve-of-growth (COG) 
analysis. For the metal-line transitions, we report rest- 
frame equivalent widths, column densities calculated us- 
ing the apparent optical depth method (AODM; Savage 
and Sembach 1991; Jenkins 1996), and statistical er- 
rors for 3<t significant transitions. At FUSE resolution, 
the majority of weak transitions are free of significant line 
saturation and the AODM approach should give accurate 
column densities. The technique has the added advantage 
that the values are free of the parametric modeling inher- 
ent to line-profile fitting. In general, upper limits reflect 
2cr statistical limits and lower limits signify probable line 
saturation. We identify each system by a redshift which 
corresponds to the centroid of the hydrogen Lyman series 
or the peak optical depth of the metal- line transitions. 



Throughout this section, the data plotted are from the 
SiC2A, LiFIA, and LiF2A channels or the STIS instru- 
ment. For the velocity plot profiles, we indicate blends 
with coincident absorption lines (primarily Galactic H2) 
as dotted orange lines. 

We caution that this section is quite detailed. The ca- 
sual reader may wish to skip to the summary table pre- 
sented at the end of the section. 




-200 -100 100 200 

Relative Velocity (km/s) 



Fig. 3. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the absorption system at z = 0.09180. In 
this and the following velocity profile plots, the dotted orange lines 
denote blends from coincident transitions and the dash-dot red line 
indicates the zero level. 



Table 5 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT z=0.09180 



Ion 


A (A) 


EW (mA) 


AODM 


-^a do pt 


HI 








14.52 ±0.05 


CIV 


1550.7700 


< 49 


< 13.70 


< 13.70 


NV 


1238.8210 


< 20 


< 13.15 


< 13.15 


OVI 


1031.9261 


73 ±8 


13.82 ± 0.05 


13.83 ± 0.04 


OVI 


1037.6167 


38 ± 10 


13.83 ±0.11 





4.1. z = 0.09180 

The lowest redshift metal-line system identified toward 
PKS0405-123 is an OVI system at z = 0.0918. Both 
transitions of the O VI doublet are detected and give a col- 
umn density iV(OVI) = 10 13 - 8±ao4 cm" 2 (Table 5). The 
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-3.4 




- 


-3.6 






-3.8 






-4.0 




/i/.-'' z=0.0918 






■V/ log N(HI) = 14.52 ± 0.04 \ 


-4.2 




b = 38±1.8km/s 



-6.5 -6.0 -5.5 

log io (ft-) 



Fig. 4. — Curve-of-growth analysis of the Lyman series for the 
gas atz = 0.09180. 



STIS+FUSE spectra reveal HI transitions for Lya, Ly/3, 
Ly7, and a marginal Lye detection (Figure 3). The results 
of a COG analysis are presented in Figure 4. Minimiz- 
ing x 2 f° r the Wa values as a function of f A, we derive 
N(tt) = io 14 - 52±0 04 cm" 2 and b = 38 ± 2 km s _1 . 

Unfortunately, the C III 977 transition for this absorber 
is blended with the Lyman series of the metal-line system 
at z = 0.167. As noted above (§ 3), CIII is particularly 
valuable for evaluating the ionization mechanism of O VI 
absorbers. With only a limited set of ionic column densi- 
ties for this absorber, both CIE and photoionization solu- 
tions are allowed although the non-detection of N V places 
tight constraints on the temperature and ionization pa- 
rameter. For collisional ionization equilibrium, the upper 
limit to A/(NV)/A(0 VI) sets a lower limit to the tem- 
perature of T > 2.5 x 10 5 K assuming [N/O] > -0.5 dex. 
This temperature is roughly consistent with the Doppler 
parameter measured for the HI gas (b ~ 40 km s _1 ). We 
also note that the O VI profile is offset by » +20 km s -1 
from the centroid of the H I profile. Perhaps the O VI gas is 
related to a more tenuous HI component. At the least, the 
kinematics raise the likelihood that this is a multi-phase 
absorber. 

For photoionization (EUVB model Q), the AT(N V) /N(0 VI) 
limit restricts the ionization parameter log t7 > [N/O] — 
1.1. Assuming [N/O] = and J 912 = 2 x 10~ 23 this im- 
plies a gas density njj < 3 x 10 cm~ 3 which is consis- 
tent with the volume density for absorbers with AT (HI) rj 
10 14 ' 5 cm -2 predicted by numerical simulations (e.g. Dave 
et al. 2001). In addition, this implies an absorber size 
£ > 30kpc. 

We can place a lower limit to the O/H metallicity by 
considering the minimum OVI/HI ionization correction. 
This gives [O/H] > — 1.4 for photoionization and [O/H] > 
—2.2 for collisional ionization. Note that these are strict 
lower limits because we have assumed all of the observed 
H I gas is associated with the O VI absorber. 

4.2. z = 0.09658 




-100 100 

Relative Velocity (km/s) 



Fig. 5. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the absorption system at z = 0.09658. 



Table 6 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT 2=0.09658 



Ion 


A (A) 


EW (mA) 


AODM 


^ adopt 


HI 








14.65 ±0.05 


CIII 


977.0200 


26 ± 7 


12.68 ±0.12 


12.68 ±0.12 


CIV 


1548.1950 


< 53 


< 13.30 


< 13.30 


NV 


1238.8210 


< 16 


< 13.06 


< 13.06 


OVI 


1031.9261 


71 ±9 


13.81 ±0.06 


13.7 ±0.2 


OVI 


1037.6167 


< 17 


< 13.63 




Silll 


1206.5000 


< 14 


< 12.02 


< 12.02 



At z = 0.09658 we identify an absorber with similar 
characteristics to the OVI system at z = 0.09180. We 
measure a nearly identical O VI column density from the 
OVI 1031 profile but note the OVI 1037 transition is 
not detected at the 3a level (Figure 5). We also caution 
that continuum placement is especially important for this 
absorber. Finally, one notes that the OVI 1031 profile 
extends to v » — 50 km s _1 whereas OVI 1037 is con- 
fined too -20kms -1 (in both the LiF2A and LiFIB 
channels). In any event, we adopt a final column density 
A(OVI) = 10 13 - 7±0 - 2 cm -2 based on the column densi- 
ties of both profiles. For now, we include this absorber in 
our sample and note that FUSE Cycle 4 observations of 
PKS0405-123 (PI: Howk) should resolve these concerns. 
A curve-of-growth analysis shows this system has nearly 
the same A(HI) value as the z — 0.09180 absorber and a 
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comparable Doppler width b = 40 ± 2 km s _1 . 

Although the C III 977 transition is blended with a co- 
incident transition, the absorption at v « +10kms _1 is 
most likely related to this absorber and we report a value 
for its column density by integrating from v ~ —10 to 
+20 km s _1 . If the C III gas arises in the same ionization 
phase as the O VI gas, the absorber cannot be explained by 
a single-phase CIE model 6 . The collisional ionization mod- 
els (Figure 2) do yield a solution with log iV(C III) /A(0 VI) 
— 1 at T w 2 x 10 5 K, yet the model also predicts 
log iV(CIV)/JV(CIII) > 1 in contradiction with the ob- 
servations. Similarly, this model would predict a NV de- 
tection unless [N/O] < —0.8 dex. In contrast, a photoion- 
ization model with logL 7 w —1.2 matches the observed 
A(CIII)/A(OVI) ratio and also predicts CIV and NV 
column densities consistent with the observed upper lim- 
its (e.g. Figure 1). Therefore, we contend photoioniza- 
tion is the principal ionization mechanism for this system. 
Adopting log £7 = —1.2 ± 0.2, we estimate the gas metal- 
licity [O/H] = -1.5 ± 0.3 dex. 




-200 -100 100 -200 -100 100 

Relative Velocity (km/s) Relative Velocity (km/s) 



Fig. 6. — Velocity profiles of the Lyman series for the partial 
Lyman limit at z = 0.16710 toward PKS0405-123. 



4.3. Z = 0.16710 

This absorption system was the principal motivation for 
the pursuit of FUSE observations of PKS0405-123. In 

6 If the absorption at v < —20 km s _1 in the O VI 1031 profile is 
confirmed in both O VI transitions (the current S/N is insufficient to 
make a definitive statement) , then one may need to consider a multi- 
phase or non-equilibrium model to explain the kinematic differences 
between the C III and O VI ions. 




« L_i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I L 

1050 1060 1070 1080 



Wavelength (A) 



Fig. 7. — Lyman limit of the absorption system at Z = 0.16710. 
The green line depicts our estimate of the flux redward of the Lyman 
limit and the blue line indicates our estimate of the flux blucward 
of the limit with a conservative error estimate (yellow lines). This 
flux decrement corresponds to log ./V(HI) = 16.45 ± 0.05. 




-8.0 -7.5 -7.0 -6.5 -6.0 -5.5 

Fig. 8. — Curve-of-growth analysis of the Lyman series for the gas 
at 2 = 0.16710. In this analysis, we have assumed two components 
with log ./V(HI) = 16.35 and 15.5 and with 40 km s _1 separation. 



particular, these data allow an analysis of the Lyman se- 
ries (Figure 6) and therefore a determination of its AT (HI) 
value. The absorber is a partial Lyman limit system, i.e., 
the optical depth at \ res t = 912A is Tgi2 < 1. Therefore, 
one can independently solve for the total A (HI) value of 
this absorber through a COG analysis and a measurement 
of T912. The A (HI) value obtained from Tgi2 is not sen- 
sitive to the Doppler parameter, instrument line-spread 
function, or any component structure within the absorber. 
In short, the T912 measurement tightly constrains the total 
HI column density with precision limited only by uncer- 
tainty in the quasar continuum. In Figure 7, we present 
the FUSE data covering the Lyman limit of the z = 0.1671 
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system. The overplotted lines mark our assessment of the 
quasar continuum redward and blueward of A = 1066A. 
Even with a very conservative estimate of the continuum 
error, the uncertainty in 7V(HI) is < 10% and we find 
iV(HI) = io 1 6- 4 5±oo5 cm -2_ 

Performing a COG analysis, we find that a single com- 
ponent model is a poor match to the observed Wa val- 
ues. Examining the low-ion profiles, we note that the C II 
and Sill profiles show two components at v\ w +5 km s _1 
and V2 ~ — 35kms _1 in Figure 7. This suggests the 
H I profile is a combination of two components and that 
a two-component COG analysis is warranted. Indeed, a 
two component solution with iVi(HI) = 10 16 35 cm -2 and 
A^HI) = 10 155 cm~ 2 and component separation Sv = 
40 km s _1 is a good description of the EW observations 
(Figure 8). In the following analysis, we will not treat the 
components separately because (i) the -ZV^HT) value is not 
well constrained and (ii) the components show similar rel- 
ative ionic ratios (with the possible exception of Nil). 

The metal-line profiles for the low and high-ion species 
from the FUSE and STIS datasets are presented in Fig- 
ures 9 and 10. The FUSE observations contribute a num- 
ber of new metal transitions, although several of the most 
valuable profiles (e.g. O II 832, O III 834, N III 989) are un- 
fortunately blended with strong Galactic transitions. The 
column density measurements and limits for all of the ob- 
served transitions are presented in Table 7. One notes 
that we report a 3er detection for the 1 988 transition 
and a > 2.5ct detection for the 01 1302 transition with 
consistent column densities. At present, we consider these 
tentative detections to be spurious primarily because the 




Relative Velocity (km/s) Relative Velocity (km/s) 



implied O metallicity is super-solar even prior to ionization 
corrections. Furthermore, the O I 1302 profile is more nar- 
row than the other low-ion profiles and is offset from the 
OI 988 profile 7 . We hope that additional, planned FUSE 
observations of PKS0405-123 (PI: Howk) will confirm (or 
contradict) this tentative OI measurement. For now, we 
adopt an upper limit to N(O ) of 10 13 85 cm~ 2 . 

Chen & Prochaska (2000) noted that no equilibrium 
photoionization or collisional ionization model could match 
the full set of ionic column densities measured for this ab- 
sorber. In particular, the ionic column densities of the 
O VI and N V ions are too large to reconcile with the low 
and intermediate ions in a single ionization phase (e.g. 
consider Figures 1,2). This conclusion is supported by the 
differences between velocity centroids and Doppler widths 
measured for the N V and O VI profiles and those of the 
other metal-line profiles. Therefore, Chen & Prochaska 
(2000) analyzed the low and intermediate ions indepen- 
dent of the O VI and N V gas and we follow this approach. 
Specifically, we examine the ionization state required to 
give the N(0 VI)/iV(N V) ratio and independently the 
physical conditions that yield the lower ions. 

Consider first the high-ions under the premise of col- 
lisional ionization. To reproduce the N(0 VI)/iV(N V) 
ratio with solar relative abundances, this implies T w 
2.5 x 10 5 K. At this temperature, this CIE phase would 
have a relatively small ionic column density for all other 
ions except CIV and SVI. The CIE model does predict 
logTV(SVI) = 12.9 which matches the column density de- 

7 The latter effect could be a small calibration error 




-200 -100 100 -200 -100 100 
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Fig. 9. — Velocity profiles of the low-ion metal transitions for the 
partial Lyman limit at z = 0.16710 toward PKS0405-123. 



Fig. 10. — Velocity profiles of the intermediate and high-ion 
metal transitions for the partial Lyman limit at z = 0.16710 to- 
ward PKS0405-123. 
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Table 7 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT z=0.16710 



Ion 


A (A) 


EW (mA) 


AODM 


N adopt 


HI 








16.45 ± 0.07 


CII 


903.6240 


125 ± 8 


> 14.24 


> 14.35 


CII 


903.9620 


158 ± 8 


> 14.13 




CII 


1334.5323 


197 ± 14 


> 14.35 




CHI 


977.0200 


374 ±8 


> 14.27 


> 14.27 


Nil 


915.6130 


104 ±7 


> 14.17 


> 14.24 


Nil 


1083.9900 


93 ± 9 


> 14.24 




NIII 


989.7990 


< 15 


< 14.59 


< 14.59 


NV 


1238.8210 


114 ± 20 


13.84 ± 0.07 


13.89 ± 0.05 


NV 


1242.8040 


78 ± 13 


13.95 ± 0.07 




01 


988.7730 


20 ± 6 


13.84 ± 0.11 


< 13.85 


01 


1302.1685 


30 ± 10 


13.74 ± 0.14 




Oil 


834.4655 


< 20 


< 13.65 


< 13.65 


OVI 


1031.9261 


377 ± 20 


> 14.74 


14.78 ± 0.07 


OVI 


1037.6167 


240 ± 23 


14.78 ± 0.07 




Sill 


1190.4158 


26 ± 7 


13.15 ± 0.09 


13.30 ± 0.04 


Sill 


1193.2897 


105 ± 11 


13.42 ± 0.05 




Sill 


1260.4221 


140 ± 11 


> 13.36 




Sill 


1304.3702 


19 ±6 


13.23 ± 0.14 




SiHI 


1206.5000 


204 ± 10 


> 13.39 


> 13.39 


SilV 


1393.7550 


140 ± 14 


13.33 ±0.04 


13.34 ±0.04 


SilV 


1402.7700 


92 ± 16 


13.44 ± 0.08 




SIII 


1190.2080 


< 13 


< 13.83 


< 13.83 


SIV 


809.6680 


< 22 


< 13.80 


< 13.80 


SIV 


1062.6620 


< 18 


< 13.89 




SVI 


933.3780 


24 ±8 


12.91 ± 0.13 


12.91 ± 0.13 


SVI 


944.5230 


< 11 


< 12.95 




Fell 


1144.9379 


< 9 


< 13.15 


< 13.15 


Felll 


1122.5260 


< 16 


< 13.60 


< 13.60 



photoionization modeling. One may formally refer to this 
metallicity as a lower limit because an unknown fraction 
of H I gas must be attributed to the more highly ionized 
gas phase. We suspect, however, that this would imply at 
most a 0.3 dex increase in [Si/H]. 



Table 8 

ELEMENTAL ABUNDANCES FOR THE ABSORBER 
AT 2=0.16710 



Ion 


[X/H] 


[X/Si+j 


C+ 


> -0.32 


> -0.03 


C++ 


> -0.88 


> -0.59 


N+ 


> 0.19 


> 0.48 


N++ 


< 0.07 


< 0.35 


N+ 4 


1.86 ±0.30 


2.15 


0° 


< 1.13 


< 1.41 


o+ 


< -1.33 


< -1.04 


O+S 


3.88 ± 0.49 


4.17 


Si+ 


-0.29 ±0.02 


0.00 


Si++ 


> -0.65 


> -0.36 


Si+ 3 


-0.04 ±0.20 


0.25 


S++ 


< 0.07 


< 0.36 


S +3 


< 0.95 


< 1.23 


S +5 


1.34 ±0.37 


1.63 


Fe+ 


< 0.78 


< 1.06 


Fe++ 


< -0.14 


< 0.15 



rived from S VI A933, although this measurement has a 
large uncertainty. Future observations of C IV AA1548, 1550 
and SVI would verify the assumption of CIE and tightly 
constrain the temperature of this phase. In passing, we 
note that the metallicity of the gas assuming T = 2.5 x 10 5 K 
and iV(Hl) = 10 16 cm- 2 is [O/H] « -ldex. 

We now focus on the lower ionization phase which we 
expect is predominantly photoionized. Figure 11 presents 
the predicted ionic column densities for a range of ions as- 
suming a plane-parallel slab of gas with iV(HI) = 10 16 5 cm -2 , 
metallicity [M/H] = -0.25 dex, and the EUVB model Q 
at z = 0.2 (see § 3). Comparisons of the model predictions 
with the observed ionic ratios constrain the ionization pa- 
rameter and provide estimates to the ionization correc- 
tions. The tightest constraints on the ionization param- — 
eter comes from the Si and N ions. The V(SiII)/iV(SiIII), 
N(N II) /N(N III) , AT(Si II) /7V(Si IV) limits and the N(Si II) /N(Si IV) 
measurement imply the dashed (lower limits), dotted (up- 
per limits) , and dash-dot (measurements) lines overplotted 
on Figure 11. Formally, no single ionization parameter sat- 
isfies all of the constraints. Given that the Cloudy calcu- 
lations are based on simplified assumptions (e.g. constant 
density, single temperature), the results are in good agree- 



Table 8 lists the [X/H] and [X/Si+] values for all of 
the ions observed in this partial LLS assuming Ubest = 
10 -3.o±o.i For iU us t r ation, we include NV and OVI to 
demonstrate these ions cannot be in photoionization equi- 
librium with the other ions because one predicts very low 
values compared to the observations. For S and Fe (marked 
by non-detections), the calculated upper limits on [S/H] 



log N(HI)=16.45! 
[M/H]=-0.25 




mcnt 8 with log U = — 2.9±0.2 dex. Applying the appropri- 
ate ionization corrections, we find [Si/H] ps — 0.25±0.2 dex 
based on the iV(SiII)/iV(HI) ratio which is relatively in- 
sensitive to U at these values. We adopt a conservative 
error on [Si/H] reflecting the systematic uncertainties of 

8 It is worth noting that adopting a softer ionizing spectrum (e.g. 
Model QG) would worsen the disagreement between the N and Si 
ions. 



Fig. 11. — Predicted ionic column densities for the absorber at 
z = 0.16710 assuming photoionization and a quasar-only extra- 
galactic UV background (model Q) at z = 0.15, a total HI col- 
umn density A r (HI) = 10 16 - 45 cm -2 , and solar relative abundances 
scaled to a metallicity [M/H] = —0.25 dex. The purple dashed line 
shows a lower limit to U based on the 7V(Si 11) /N (Si III) upper limit 
whereas the dotted purple lines show upper limits to U based on 
limits to the iV(NII)/iV(NHI) and jV(SiII)/AT(SiIV) ratios. Fi- 
nally the N (Si II)/ N (Si IV) value and la error constrain U to the 
region indicated by the dash-dot purple lines. 
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and [Fe/H] are consistent with expectation. The C ions 
are also in reasonable agreement although the lower limits 
from C II nearly conflict with the Si abundance under the 
presumption of solar relative abundances. As noted above, 
if the O I transitions are true detections they imply a large 
O abundance and a correspondingly large O/Si ratio. If 
future observations yield an N(0°) value near our upper 
limit, the super-solar O/Si ratios might indicate gas en- 
riched by massive (M > 15M Q ) Type II SN (Woosley & 
Weaver 1995). 

Because the FUSE+STIS datasets provide full coverage 
of the dominant Si and N ions indicated by the photoion- 
ization modeling, it is possible to calculate the N/Si ratio 
largely independent of ionization corrections for the pho- 
toionized gas 9 : 

N _ TV (Nil) + V(NHT) 

Si ~ JV(Sill) + iV(Silll) + iV(SilV) ' (1) 

Summing over the column densities and assuming reason- 
able values for the limits, we find log(N/Si) > +0.5 or 
[N/Si] > +0.2 dex. We derive a similar value by consid- 
ering only the iV(NII)/./V(SiII) ratio which has a small 
dependence on U for logf < —2.7 dex. This N/Si over- 
abundance is a surprising result; there are few examples 
of enhanced N abundances in the local universe. We also 
note that the O and C observations are consistent with 
enhancements of these elements relative to Si suggesting 
an overabundance of all the light elements. One of the few 
sites where large N abundances are observed are quasars 
(Osmer 1980; Bentz, Hall, & Osmer 2004). Given the 
properties of high mctallicity, enhanced light elements, and 
substantial separation from any bright galaxy (Prochaska 
et al. 2004), one may speculate that this gas was related 
to AGN activity. 

4.4. z = 0.18250 and z = 0.18290 

The system at z = 0.183 is notable for showing large 
N(0 VI)/iV(HI) variations in a pair of absorbers sepa- 
rated by 6v < 100 km s _1 . Figure 12 plots the HI and 
metal-line profiles for the two absorbers with z — 0.18250 
(v = -25 km s- 1 ) and z = 0.18290 (v = +50 km s" 1 ). We 
performed separate COG analyses and derived H I column 
densities iV(HI) = 10 14 - 50±0 05 cm" 2 and 10 14 08±01 cm" 2 
and Doppler parameters b = 27 ± 1 km s -1 and b = 26 ± 
5 km s _1 respectively. The only metal-line absorption that 
we confidently identify with this system is O VI gas related 
to the z — 0.1829 absorber. There is a hint of OVI ab- 
sorption at the velocity corresponding to the z — 0.1825 
absorber, but the measured equivalent width is less than 
a 3a detection. At this redshift, the FUSE data provide 
coverage of the O III 832 and O IV 787 transitions of this 
absorber. Unfortunately, the O IV 787 transition is signifi- 
cantly blended with an H2 transition. The O III 832 region 
is also partially contaminated but provides a meaningful 
upper limit on the OIII column density (Tables 9, 10). 

We detect no other ions for these two absorbers. Most 
importantly, we place an upper limit to N(C III) < 10 12 ' 4 cm 
for the absorber at z = 0.18290 10 . If the gas is photoion- 

9 Note jV(NlV) and AT(NV) are assumed to be smaller than the 
upper limit obtained for jV(N III) 

10 The CIII 977 profile is blended with an H2 transition at the 
velocity corresponding to 2 = 0.18250. 



Table 9 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT ^=0.18250 



Ion 


A (A) 


EW (mA) 


AODM 


■^V 0, do pt 


HI 








14.90 ±0.05 


NV 


1238.8210 


< 14 


< 12.99 


< 12.99 


on 


834.4655 


< 21 


< 14.55 


< 14.55 


OIII 


832.9270 


< 18 


< 13.73 


< 13.73 


OIV 


787.7110 


< 19 


< 13.83 


< 13.83 


OVI 


1031.9261 


< 28 


< 13.77 


< 13.77 


OVI 


1037.6167 


< 23 


< 13.95 




Sill 


1260.4221 


< 15 


< 12.21 


< 12.21 


Silll 


1206.5000 


< 14 


< 11.99 


< 11.99 


SilV 


1393.7550 


< 28 


< 12.64 


< 12.64 
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IONIC COLUMN DENSITIES FOR THE ABSORBER 






AT 2=0.18290 




Ion 


A (A) 


EW (mA) 


AODM 


N adopt 


HI 








14.08 ±0.10 


CIII 


977.0200 


< 11 


< 12.42 


< 12.42 


NV 


1238.8210 


< 15 


< 13.05 


< 13.05 


OIII 


832.9270 


< 25 


< 13.79 


< 13.79 


OIV 


787.7110 


< 22 


< 14.61 


< 14.61 


OVI 


1031.9261 


72 ± 18 


13.94 ±0.12 


13.96 ±0.07 


OVI 


1037.6167 


44 ±9 


13.97 ±0.09 




Silll 


1206.5000 


< 14 


< 12.01 


< 12.01 


SilV 


1393.7550 


< 27 


< 12.68 


< 12.68 




CIII 977 



-100 100 

Relative Velocity (km/s) 



Fig. 12. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the pair of absorption systems at 2 = 0.183. 
The dotted line in the figure corresponds to 2 = 0.18270. 
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ized, the N(0 VI)/iV(CIII) ratio requires an ionization 
parameter log U > — 1 for solar relative abundances which 
implies [O/H] > — 1 (e.g. Figure 1). Under the assumption 
of collisional ionization, however, the iV(CIII)/iV(0 VI) 
limit requires T > 2 x 10 5 K and suggests metallicity 
—2 < [O/H] < — 1. If we ignore the Doppler parame- 
ter of the HI gas which is uncertain, both solutions are 
allowed but imply temperatures for the gas which differs 
by more than one order of magnitude. A detection or 
sensitive limit to the iV(CIV) value would likely resolve 
the degeneracy; this measurement could be achieved with 
future HST+STIS observations of PKS0405-123. 

Irrespective of the ionization mechanism, the absorbers 
are notable for exhibiting N (O VI) / N (HI) ratios differ- 
ing by more than an order of magnitude in gas sepa- 
rated by Sv < 100 km s" 1 (see also Savage et al. 2002). 
Under the assumption of photoionization, the variation 
in N(0 VT)/iV(HI) can only be explained through large 
differences in [M/H] and/or nn- For example, the z = 
0.18250 absorber may have substantially higher gas den- 
sity and a correspondingly lower ionization parameter. In 
this case, however, the models predict iV(OIV) values in- 
consistent with the observed upper limits. Therefore, if 
the gas is photoionized the absorber at z = 0.18290 has 
at least 10 x higher metallicity. This would indicate a re- 
markable level of inhomogeneity in metal enrichment and 
would be particularly surprising given the z = 0.18290 ab- 
sorber has lower HI column density. If collisional ioniza- 
tion is the dominant mechanism, then the N (O VI) / N (HI) 
ratios imply large temperature differences if the two ab- 
sorbers have comparable metallicity (e.g. T > 10 6 K for 
the gas at z = 0.18250 and T w 3 x 10 5 K for the gas 
at z = 0.18290). If we require T > 2 x 10 5 K for the 
z = 0.18290 system based on the non-detection of CIII, 
then we derive a temperature for the z = 0.18250 absorber 
which is inconsistent at greater than 5<r significance from 
the COG analysis. Therefore, it is likely that a CIE model 
would also require a large (> 0.3 dex) difference in [O/H]. 

4.5. z = 0.36080 

There is a complex, highly ionized absorption system 
at z = 0.36080 toward PKS0405-123. The FUSE+STIS 
datasets show Lya-Ly6 absorption (Figure 13) which shows 
a hint of two closely separated components. A single com- 
ponent COG analysis, however, provides a good descrip- 
tion of the data and reveals AT (HI) = 10 1512±0 05 cm" 2 
and b — 27 ± 1 kms -1 . We identify a series of metal- 
line transitions including absorption from O III O IV, and 
CIII (Figure 14, Table 11). We also detect NIV 765 and 
Si III 1206 but report upper limits because of line blending. 
Finally, there is a 2a detection of the N III 989 transition 
and a feature at the expected position of NV 1238. We 
suspect the latter feature is related to fixed pattern noise or 
detector artifact because the observed optical depth is only 
marginally consistent with the null detection of N V 1242 
and because of the implied N abundance: [N/Si] > 0.5 dex. 
Interestingly, this absorber shows no low-ion gas and also 
no OVI absorption. We predict the absorber will show 
significant C IV absorption and, therefore, may be an ana- 
log to high z C IV systems with comparable H I column 
density (e.g. Sargent, Steidel, & Boksenberg 1988). In- 
deed, this is apparently the case as Williger et al. (2004) 
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Fig. 13. — Velocity profiles of the Lyman series for the absorption 
system at z = 0.36080 toward PKS0405-123. 
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Fig. 14. — Velocity profiles of the metal-line transitions for the 
absorption system at z = 0.36080. 
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find log iV(CIV) = 13.76 ± 0.16 based on their analysis of 
STIS G230M archival data. 



Table 11 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT z=0.36080 



Table 12 

ELEMENTAL ABUNDANCES FOR THE ABSORBER 
AT z=0.36080 



Ion 


A (A) 


EW (mA) 


AODM 


$ adopt 


HI 








15.12 ±0.05 


CII 


1036.3367 


< 13 


< 13.21 


< 13.21 


cm 


977.0200 


152 ±7 


> 13.78 


> 13.78 


Nil 


1083.9900 


< 19 


< 13.47 


< 13.47 


NIII 


989.7990 


< 11 


< 13.31 


< 13.31 


NIV 


765.1480 


< 8 


< 13.22 


< 13.22 


NV 


1238.8210 


< 23 


< 13.52 


< 13.52 


NV 


1242.8040 


< 32 


< 13.75 




on 


834.4655 


< 10 


< 13.25 


< 13.25 


OIII 


832.9270 


72 ±5 


> 14.17 


> 14.17 


OIV 


787.7110 


75 ±7 


> 14.26 


> 14.26 


OVI 


1031.9261 


< 16 


< 13.30 


< 13.30 


Sill 


1260.4221 


< 28 


< 12.55 


< 12.55 


SiHI 


1206.5000 


< 19 


< 12.64 


< 12.64 


Felll 


1122.5260 


< 19 


< 13.67 


< 13.67 



Examining the collisional ionization solutions described 
by Figure 2, a solution with T f» 1.3 x 10 5 K matches 
the relative column densities of the O ions. We do not fa- 
vor this scenario, however, for several reasons. First, col- 
lisional ionization implies iV(C UT)/N(0 III) w -0.8 yet 
we observe 7V(CIII)/7V(OIII) > —0.4 under the assump- 
tion that the C III 977 profile requires larger corrections 
for line-saturation than OIII 832. Second, the tempera- 
ture is inconsistent at > 5er with the value implied by the 
measured H I Doppler parameter. Third, the predicted 
AT(NIII)/iV(NlV) ratio is only marginally consistent with 
the observed limits. Although none of these arguments is 
decisive, together they imply collisional ionization is not 
the dominant ionization process. A conclusive test (e.g. 
observations of the CIV doublet) would be valuable as 
CIE models imply significantly lower metallicity and cor- 
respondingly higher iV(H) than photoionization. 

Assuming photoionization, the measurements and lim- 
its on the column densities of the O ions tightly con- 
strain the ionization state of this gas (Figure 15). Treat- 
ing the OIII and OIV column densities as lower limits 
because of line saturation, the resulting N(0 II)/iV(0 IV) 
and N(0 III) /N(0 VI) limits require -2.1 < log U < -1.7. 
If we adopt log[JV(0 m)/N(0 IV)] = 0.1 ±0.2, this places 
a similar constraint on the U parameter. In Table 12, 
we list the ions analyzed for this absorber and [X/H] val- 
ues corresponding to Ubest — 10 — 1 ° -2 . The principal 
result is that the lower limits to O III and O IV imply 
[O/H] > —0.7, an enrichment level consistent with the 
lower limit to the C III column density. 

Because the ionization corrections for OIII/CIII are 
nearly constant at — 0.2dex for the relevant ionization pa- 
rameter, we estimate [O/C] < +0.2 based on the satu- 
rated OIII A832 and CIII A977 profiles. In this absorber, 
at least, it is unlikely the gas has a super-solar O/C ra- 
tio, in contrast to the abundances inferred for a sample of 
z > 1 absorption systems with comparable ionization state 
(Vogel & Reimers 1995). The upper limit measurements 
for the nitrogen ions also place a limit on the N/O abun- 
dance. In particular, note that the predicted NIII column 



Ion 


rv /m 

[A/ri] 


[v + l 


r<+ 


<. U.Uo 


V U. to 


C++ 


> —1.05 


> —0.30 


N+ 


< 1.08 


< 1.83 


N++ 


< -0.96 


< -0.21 


N+3 


< -0.92 


< -0.17 


N+4 


< -0.07 


< 0.68 


o+ 


< 0.15 


< 0.90 


o++ 


> -0.75 


> 0.00 


0+3 


> -0.76 


> -0.01 


Q+5 


< -0.54 


< 0.20 


Si+ 


< 1.34 


< 2.09 


Si++ 


< 0.12 


< 0.87 


Fe++ 


< 3.77 


< 4.52 



density closely traces A^OIII) for logU « —2. Therefore, 
the ratio of these ions provides a measurement of N/O 
largely independent of uncertainties in the ionization pa- 
rameter: [N/O] < -0.2 dex. Noting that the AT(NIII) 
value is likely to be within rs 0.1 dex of the adopted upper 
limit (weak absorption is apparent at the expected posi- 
tion of NIII 989), we report [N/O] = -0.3 ± 0.15 dex. 
This relative abundance follows the [N/O] vs. [O/H] trend 
observed for galaxies in the local universe (e.g. Henry, Ed- 
munds, & Koppen 2000). 

Because this absorber is characterized by higher ion- 
ization states, the results are sensitive to the assumed 
shape of the EUVB radiation field. If we adopt the softer 
model QG, the observed ionic ratios for oxygen require a 
higher ionization parameter (U w —1.4) and imply lower 



N(HI)=15.10 
[M/H]=-0.50 
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Fig. 15. — Predicted ionic column densities for the absorber at z = 
0.36080 assuming photoionization and a quasar-only extragalactic 
UV background (model Q) at z = 0.35, a total HI column density 
Af(HI) = 10 15,10 cm -2 , and solar relative abundances scaled to a 
metallicity [M/H] = —0.50 dex. The purple dashed line shows a 
lower limit to U based on the N(0 ll)/N(0 IV) upper limit whereas 
the dotted purple line shows an upper limits to U based on the lower 
limit to N(0 IV) /N(0 VI). Finally the N(0 III) /N(0 IV) value and 
1<T error constrain U to the region indicated by the dash-dot purple 
lines. 
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elemental abundances ([O/H] > —1.2). The latter result 
is a consequence of the higher H ionization fraction of the 
softer spectrum. Again, a measurement of the CIV col- 
umn density would help resolve this ambiguity. Addition- 
ally, tighter constraints on the column densities of OVI 
and the N ions would probe the shape of the ionizing spec- 
trum. 



Table 13 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT z=0.36332 



Ion 


A (A) 


EW (mA) 


AODM 


^ adopt 


HI 








13.43 ±0.10 


CII 


1036.3367 


< 13 


< 13.22 


< 13.22 


cm 


977.0200 


22 ± 7 


12.64 ±0.12 


12.64 ±0.12 


Nil 


1083.9900 


< 15 


< 13.33 


< 13.33 


NIII 


989.7990 


< 19 


< 13.47 


< 13.47 


NIV 


765.1480 


42 ±7 


> 13.22 


> 13.22 


NV 


1238.8210 


< 27 


< 13.34 


< 13.35 


NV 


1242.8040 


< 31 


< 13.70 




on 


834.4655 


< 13 


< 13.33 


< 13.33 


OIII 


832.9270 


< 11 


< 13.40 


< 13.40 


OIV 


787.7110 


28 ±6 


13.73 ±0.10 


13.73 ±0.10 


OVI 


1031.9261 


23 ± 8 


13.44 ±0.11 


13.44 ±0.11 


OVI 


1037.6167 


< 11 


< 13.45 




NeVIII 


770.4090 


< 13 


< 13.55 


< 13.55 


Sill 


1260.4221 


< 31 


< 12.53 


< 12.53 


Silll 


1206.5000 


< 24 


< 12.26 


< 12.26 



4.6. z = 0.36332 
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There is a neighboring absorber (at +550 km s 1 ) to 
the z = 0.36080 system which shows OIV, OVI, NIV, 
and C III absorption yet a weak and unusual Lya profile 
(Figure 16). In particular, one notes that the metal-lines 
are significantly offset from the H I absorption. Integrating 
the Lya profile from v = —30 to +80 km s -1 , we derive 
an HI column density AT (HI) = 10 13 4±a2 cm~ 2 . In the 
following, we will adopt this value but one may consider it 
an upper limit to the HI column density associated with 
this metal-line system. It would be extremely valuable 
to have higher S /N of the C IV and Lya spectral regions 
for this absorber to more confidently compare their line 
profiles. We suggest this system is an example of an O VI 
absorber without detected neutral hydrogen gas which im- 
plies a very high temperature and/or metallicity. Table 13 
lists the ionic column densities for this absorber which in- 
cludes an upper limit on A (NIV) based on the blended 
NIV 765 profile. 



Table 14 

ELEMENTAL ABUNDANCES FOR THE ABSORBER 
AT 2=0.36332 



Ion 


[X/H] 


[X/Si+] 


C+ 


< 2.36 


< 2.34 


C++ 


-0.37 ±0.10 


-0.39 


N+ 


< 3.22 


< 3.20 


N++ 


< 0.97 


< 0.95 


N+3 


> 0.34 


> 0.32 


N+ 4 


< 0.58 


< 0.56 


o+ 


< 2.59 


< 2.57 


o++ 


< 0.31 


< 0.29 


0+3 


0.02 ±0.07 


0.00 


o+ 5 


0.04 ± 0.44 


0.02 


Ne+ 7 


< 0.05 


< 0.03 


Si+ 


< 4.79 


< 4.77 


Si++ 


< 2.92 


< 2.90 



Fig. 16. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the absorption system at z = 0.36332. 



It is surprising, however, that the absorber also shows 
C III 977. Regarding collisional ionization equilibrium mod- 
els, there is no single temperature solution which would 
give A(CIII)/A(OVI) « 0.15 (T < 2 x 10 5 K) and 
A(OIV)/V(OVI) w 2 (T « 2.5 x 10 5 K) assuming solar 
relative abundances. Although the difference in these tem- 
peratures may appear small, the relative abundances of 
A(C III) /A+O VI) and A(OIV)/A(0 VI) are very sensi- 
tive to temperature (see also Howk et al. 2004) . For exam- 
ple, assuming T = 2.5 x 10 5 K would underpredict [C/O] 
by two orders of magnitude. Therefore, a single-phase 
CIE model cannot fully describe this absorber. Adopt- 
ing photoionization as the dominant mechanism, the best 
model has log U = —1.4+0.2 assuming the EUVB model Q 
and [M/H] =0. Table 14 presents the [X/H] and [X/0+ 3 ] 
values for U best = 10" 14±0 ' 2 . Even adopting A(HI) = 
10 13 4 cm -2 , we calculate a solar oxygen abundance for 
this absorber. Independent of the 7V(HI) value we find 
[O/C] w +0.4 + 0.2 dex. These abundances are compara- 
ble to those derived from the nearby z = 0.36080 system. 
Although the velocity separation is too large to associate 
spatially the gas, the 'clouds' may have a common enrich- 
ment history suggesting a similar physical origin. 
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Table 15 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT z=0.49510 



Ion 


A (A) 


EW (mA) 


AODM 


N adopt 


HI 








14.39 ± 0.07 


CII 


903.9616 


< 20 


< 13.08 


< 13.08 


cm 


977.0200 


61 ± 7 


13.18 ± 0.06 


13.18 ± 0.06 


Nil 


1083.9900 


< 22 


< 13.47 


< 13.47 


NIII 


989.7990 


< 12 


< 13.29 


< 13.29 


NIV 


765.1480 


< 16 


< 12.88 


< 12.88 


on 


834.4655 


< 13 


< 13.37 


< 13.37 


OIII 


832.9270 


< 15 


< 13.59 


< 13.59 


OIV 


787.7110 


105 ±9 


> 14.34 


> 14.34 


ov 


629.7300 


164 ±9 


> 14.41 


> 14.41 


OVI 


1031.9261 


133 ±9 


14.25 ±0.04 


14.26 ±0.03 


OVI 


1037.6167 


94 ± 10 


14.31 ±0.06 




NeVIII 


770.4090 


26 ± 7 


13.75 ±0.10 


13.75 ±0.10 


NeVIII 


780.3240 


< 14 


< 13.86 




MgX 


624.9500 


< 16 


< 14.22 


< 14.22 


SVI 


933.3780 


< 9 


< 12.62 


< 12.62 



4.7. z = 0.4951 

The highest redshift metal-line system along the sight- 
line to PKS0405-123 is at z = 0.495 and exhibits detec- 
tions of GUI, OIV, O V, and possibly NeVIII transitions. 
Williger et al. (2004) performed a profile fit to the Ly/3 
transition and measured logiV(HI) = 14.39 ± 0.07 and 
b = 60±15 km s _1 consistent with our W\ measurements. 
The HI and metal-line profiles are presented in Figure 17 
and the ionic column densities are listed in Table 15. 

A full analysis of this absorber is presented in Howk 



et al. (2004) and we only summarize a few key results 
here. The FUSE spectra provide constraints on a num- 
ber of transitions which are very rarely observed in inter- 
vening quasar absorption line systems including OV 629, 
NIV 765, and the NeVIII and MgX doublets. Altogether, 
the FUSE±STIS datasets provide a comprehensive analy- 
sis of the ionization state of this gas. While a collisional 
ionization model with T « 2.3 x 10 5 K can reproduce the 
observed O ionic column densities, this model underpre- 
dicts AT(CIII) by over 2 orders of magnitude. Therefore, 
Howk et al. (2004) rule out a single-phase CIE model. In 
contrast, the gas is well modeled by a single-phase pho- 
toionization model with \ogU « —1.2 although a two- 
phase model (collisional plus photoionization) is also al- 
lowed. We summarize the physical properties of this gas 
in the following section. 

5. N (HI) > 10 14 cm" 2 ABSORBERS WITHOUT METAL-LINE 
TRANSITIONS 

In our analysis of PKS0405-123 we have searched for 
metal-line transitions in every Lya absorption system iden- 
tified by Williger et al. (2004), with special attention to 
AT (HI) > 10 14 cm~ 2 systems. This allows us to investi- 
gate the metallicity distribution of a 'complete' sample of 
low redshift, AT(HI) > 10 14 cm -2 Lya clouds. Perhaps the 
most remarkable example is the absorber at z = 0.40570. 
Although this absorber shows nearly the same HI col- 
umn density as the metal-line system at z = 0.36080, 
it is remarkable for exhibiting no metal-line transitions. 
In Figure 18 we present the Lyman series and a set of 
the strongest metal transitions covered by the FUSE and 
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Fig. 17. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the absorption system at z = 0.49510. 



Fig. 18. — Velocity profiles of the Lyman series and metal-line 
transitions analyzed for the absorption system at z = 0.40570. Note 
that none of the metal-line transitions show significant detections. 
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STIS spectra (Table 16). A COG analysis gives 7V(HI) = 
iQi4.85±o.o5 cm -2 with Doppler parameter b = 34±2 km s _1 . 
Formally, there is a detection of OVI gas via the broad 
O VI 1031 transition, yet the continuum varies significantly 
near 1450A in the STIS spectrum and the measured value 
is inconsistent with the upper limit implied by the O VI 1037 
transition. In the following, we adopt the upper limit from 
the O VI 1037 transition. 



Table 16 

IONIC COLUMN DENSITIES FOR THE ABSORBER 
AT 2=0.40570 



Ion 


A (A) 


EW (mA) 


AODM 


N adopt 


HI 








14.85 ±0.05 


CII 


903.9616 


< 13 


< 12.92 


< 12.92 


CII 


1036.3367 


< 13 


< 13.23 




cm 


977.0200 


< 14 


< 12.53 


< 12.53 


Nil 


1083.9900 


< 16 


< 13.32 


< 13.32 


NIII 


989.7990 


< 16 


< 13.42 


< 13.42 


on 


834.4655 


< 14 


< 13.37 


< 13.37 


OIII 


832.9270 


< 15 


< 13.50 


< 13.50 


OIV 


787.7110 


< 12 


< 13.43 


< 13.43 


OVI 


1031.9261 


< 14 


< 13.63 


< 13.63 


OVI 


1037.6167 


< 19 


< 13.67 




SiHI 


1206.5000 


< 32 


< 12.34 


< 12.34 



Assuming an equilibrium photoionization model with 
EUVB model Q, the upper limits to the observed ions 
constrain the metallicity of the g function of the ion- 

ization parameter. In Figure 19, we plot the upper limit 
to the gas metallicity [M/H] for a range of U parameters. 
The plot symbols indicate the ion which places the most 
sensitive limit to [M/H]. For \ogU w —3, the gas must 
have [M/H] < —1, but this would imply rin ~ 10~ 3 cm~ 3 
(for Jg 12 = 2 x 10~ 23 ) which is considerably higher than 
expected for such a low A(HI) absorber. A more realistic 
value is logC > -2 for iV(HI) = 10 14 8 cm" 2 (§ 6.3) which 
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Fig. 19. — Upper limit to the metallicity of the absorber at 
z = 0.40570 as a function of ionization parameter. At each U value, 
the point style indicates the ion with the tightest constraint on the 
metallicity. For most of the parameter space, C III places the low- 
est limit to [M/H]. The dashed line indicates the likely ionization 
parameter for an absorber with JV(HI) = 10 14 8 cm -2 (§ 6.3). 



implies [M/H] < — 2dex. The results indicate that this 
absorber has 5 to 100 times lower metal enrichment than 
the otherwise similar absorption system at z = 0.36080. 
Unless this absorber is a special example, systems with 
AT(HI) « 10 15 cm~ 2 have large metallicity dispersion at 
these redshifts. We also note that a CIE model would 
imply even lower metallicities then the results shown in 
Figure 19. 

As an aside, we stress that C III places the tightest limit 
on the metal abundance for the largest range of parameter 
space. This ion and its corresponding C III 977 transition 
deserve greater attention in discussions of the Lya for- 
est. At higher redshift, the transition is frequently blended 
within the forest, yet we emphasize that it is a very pow- 
erful diagnostic of metal enrichment and ionization con- 
ditions within absorbers like the one at z = 0.40570 (c.f. 
Schaye et al. 2003). 

In addition to the absorber at z = 0.4057, there are sev- 
eral absorbers with no unambiguous, statistically signifi- 
cant metal-line transitions. These are: z = 0.030, AT (HI) = 
10 14 - 2 cm" 2 ; z = 0.351, AT (HI) = 10 142 cm" 2 ; 
z = 0.409, AT(HI) = 10 14 - 4 cm- 2 ; z = 0.538, A(HI) = 
10 14 ' 2 cm -2 . In several cases, there are suggestions of 
metal-line absorption but line blending or poor S/N pre- 
vents a definitive measurement. In all cases, one can es- 
timate an upper limit to their metal enrichment with an 
analysis similar to the one we performed for the z = 0.40570 
system. Under the assumption of photoionization and as- 
suming logf/ > —2, we find [M/H] < — 1.5 based on the 
non-detection of CIII 977 and/or OVI and adopt this up- 
per limit in the following discussion. 

6. DISCUSSION 

Section 4 presented a detailed analysis of the metal- 
line systems observed along the sightline to PKS0405-123. 
With the combined FUSE+STIS datasets, we have ob- 
tained constraints on the ionization mechanism, temper- 
ature, metallicity, and elemental abundance ratios of the 
gas. We summarize these properties in Table 17 which 
lists the redshift, A(HI) value, N(0 VI) value, the ioniza- 
tion parameter, metallicity [M/H] p ^ ot , the total H column 
density A(H) assuming photoionization, the mean volume 
density Uh, and an estimate of the length scale I assum- 
ing photoionization and J912 = 2 x 10~ 23 cgs. The table 
also lists the temperature and metallicity [M/H] co j; for col- 
lisionally ionized gas, and an estimate of N/a, the abun- 
dance of nitrogen relative to an a-element (e.g. Si, O). The 
following sub-sections consider various implications of our 
analysis on the nature of metal-line systems in the low 
redshift universe. 

6.1. Tracing the WHIM with O VI Absorbers 

As noted in the introduction, current theoretical expec- 
tation is that a significant fraction of baryons at z < 1 are 
in a warm-hot intergalactic medium (WHIM). One of the 
few means of detecting this gas is through O VI absorp- 
tion. To this end, Tripp and collaborators have recently 
surveyed OVI gas along a number of quasar and AGN 
sightlines (Tripp et al. 2000; Savage et al. 2002). Their 
latest results indicate an incidence dNovi/dz = 14+g for 
an Wa limit of 50mA. This implies a baryonic mass den- 
sity f?ovi > 0-002hj^ assuming the average metallicity 
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Table 17 
SUMMARY TABLE 



z 


iV(HI) 
(cm z ) 


bni 


iV(OVI) 
( cm ) 


Ion a 


logLV 


[M/R\ phot 


log N(H) 
(cm z ) 


(lO" 5 ) 


l c 
(kpc) 


T co u 
(10 5 K) 


[M/Hj co/; [N/aj 


0.09180 


14.5 


38 


13.8 


?? 


> -1.5 


> -1.4 


18.8 


< 2 


> 90 


> 2.5 


> -2.2 


0.09658 


14.7 


40 


13.7 


Photo 


-1.2 


-1.5 


19.2 


1.2 


440 






0.16710 


16.5 


38 


14.8 


Multi 


-2.9 


-0.25 


18.8 


60 


3 




>0.2 


0.18250 


14.9 


27 


< 13.8 


?? 










< -1 






0.18290 


14.1 


26 b 


14.0 


Coll 


> -1 


> -1 


> 18.4 


3 


« -1.5 






0.36080 


15.1 


27 


< 13.3 


Photo 


-2.0 


> -0.7 


18.6 


8 


16 




-0.3 


0.36332 


13.4 




13.4 


Photo 


-1.4 





17.4 


2 


1 






0.49510 


14.4 


60 b 


14.3 


Photo 


-1.3 


> -0.3 


18.5 


1.7 


60 


2.6 


> -1 < -O.f 



"This column lists the dominant ionization mechanism for the gas as inferred from the metal-line ratios. In several cases, the situation is 
ambiguous or multiple processes are required. 
6 Therc is a > 25% uncertainty in this value. 
c These values assume J912 = 2 X 10~ 23 cgs. 



is 0.1 solar and a very conservative ionization correction 
{N(OVl)/N(0) < 0.2). 

With the large redshift of PKS0405-123 and its exten- 
sive UV spectroscopy, we can derive OVI statistics for 
this sightline. We have detected six O VI systems to an 
W\ limit of 30 mA over a non-contiguous redshift path of 
Az — 0.38. This pathlength was determined by identifying 
all of the regions free of Galactic H 2 and other coincident 
transitions with sufficient S/N to identify and measure 
both members 11 of the OVI doublet to a 3cr Wa thresh- 
old of 60 mA. We also searched for any O VI absorbers 
without corresponding Lya absorption (to the same limit) 
and found no examples. Assuming Poisson statistics (e.g. 
Gehrels 1986), we find dN vi/dz = 16±l (±la) which is 
consistent with Savage et al. (2002). If we combine our 
results with Savage et al. (2002) - ignoring the small dif- 
ference in sensitivity - we find: dNovi/dz = 15^ (±ler). 

An important aspect of the present work is that our ob- 
servations place constraints on the ionization state and/or 
metallicity of several OVI absorbers. For the systems 
at z = 0.09650, z = 0.36332 and z = 0.49510, the de- 
tection of C III 977 precludes a single-phase CIE model 
and therefore suggests the gas is photoionized, in multiple 
ionization phases, or in a non-equilibrium state (see also 
Howk et al. 2004). If these absorbers have low temper- 
ature (T < 10 5 K), they are unlikely to contribute to the 
WHIM or present a significant reservoir of baryons. Fu- 
ture observations of the C IV profiles of these absorbers 
would be particularly valuable to address this issue. In 
contrast, the OVI gas in the absorbers at z = 0.16710 and 
z = 0.18290 are well modeled by a single-phase CIE model 
and we believe this gas is collisionally ionized. Finally, the 
absorber at z = 0.09180 is well modeled by both CIE and 
photoionization solutions. If the fraction of photoionized 
O VI absorbers is roughly half, then flo vi in the WHIM 
is roughly 50% the value that one would infer assuming 
all of the O VI gas is collisionally ionized. It will be im- 

11 Note that these selection criteria exclude the likely OVI sys- 
tem at z = 0.08 because the OVI 1037 transition is blended with a 
coincident transition. 



portant to examine the environment of the absorbers as a 
function of ionization state, e.g. to determine if the colli- 
sionally ionized gas arises preferentially in galactic halos 
(Scmbach et al. 2003) or large-scale structures. 

Our analysis also constrains the metallicity of the gas 
giving rise to O VI absorption. For the absorbers pre- 
sumed to be photoionized, the gas metallicity ranges from 
[M/H] « —1.5 to dcx. The collisionally ionized gas shows 
systematically lower metallicity ([O/H] « —2 to — ldex) 
because this gas has much larger hydrogen ionization frac- 
tion. Altogether, these measurements are consistent with 
the values assumed and expected for O VI systems if the 
gas is related to the WHIM (e.g. Savage et al. 2002; Dave 
et al. 2001). It is worth noting, however, that several of 
the OVI lines are at the detection limit of this dataset 
and, therefore, absorbers with significantly lower metallic- 
ity would not have been detected. Unfortunately, it may 
not be possible with FUSE or HST/STIS to probe metal- 
licities below [O/H] w -2. 

6.2. Tracing the WHIM with Ne VIII Absorbers 

At temperatures T > 5 x 10 5 K, O VI is no longer the 
dominant ion of oxygen. This is unfortunate because cur- 
rent simulations of the WHIM predict the bulk of the gas 
is at temperatures T > 10 6 K. For these reasons, X-ray 
spectroscopy of O VII and O VIII transitions offer a more 
direct probe of the WHIM (e.g. Fang et al. 2002) yet 
current technology allows a search for OVII and OVIII 
in only the brightest few X-ray sources. Even with the 
next generation of X-ray telescopes, surveys for O VII and 
OVIII will have limited impact. An alternative means 
of probing collisionally ionized gas with T > 10 6 K is 
through the NeVIII AA770,780 doublet. At these tem- 
peratures, NeVIII is the dominant Ne ion and this rel- 
atively strong pair of transitions can be observed with 
FUSE and HST/STIS spectroscopy at z > 0.2. Specifi- 
cally, FUSE spectroscopy with comparable spectral qual- 
ity to the observations presented here will be sensitive to 
NeVIII gas in absorbers with logiV(HI) + [M/H] > 13 for 
10 5 5 K < T < 10 7 K. 
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For sightlines with significant Galactic H2 absorption 
(PKS0405-123 is a moderate example), a survey for Ne VIII 
is difficult for A < 1000 A or z < 0.3. Nevertheless, we 
have conducted a search for Ne VIII absorption by exam- 
ining the NeVIII regions for all Lya absorbers detected 
along the PKS0405-123 sightline. We found no signifi- 
cant detections. Furthermore, for A > 1000 A we have 
found no pairs of absorption features with the appropri- 
ate separation of the NeVIII doublet, independent of the 
presence of a corresponding Lya feature. Altogether, we 
estimate the non-contiguous pathlength searched to be 
AzncVUI = 0.11. Assuming Poisson statistics, we set a 
conservative 95% upper limit to the incidence of NeVIII 
absorption: dN$ievm/dz < 40 for a 3a rest-EW limit of 
» 40mA. This result will be extended to lower Wa limits 
and pursued along additional sightlines with forthcoming, 
approved FUSE programs (PI: Howk). We will also ex- 
tend the analysis to include a search for Mg X absorption 
for systems with z > 0.5. 

6.3. Ionization of the IGM: A Possible A(H) Conspiracy 

The analysis described in § 4 allows us to investigate 
global trends in the ionization state of low z metal-line 
systems. One key result is that a significant fraction of 
the systems are predominantly photoionized. This is not 
a surprising result; current expectation is that the ma- 
jority of the low z Lya forest is photoionized gas (e.g. 
Dave & Tripp 2001). Unlike studies of the majority of 
Lya 'clouds', however, an analysis of the metal-line sys- 
tems provides an assessment of the ionization fraction and 
metallicity of this gas. 

Focusing on the subset of systems where observations 
constrain the ionization state, we can examine trends be- 
tween U and A (HI). Figure 20 presents a plot of logU vs 
log A (HI) for all of the systems where we believe photoion- 
ization is the dominant mechanism. We caution, again, 
that in several cases the presumption of photoionization is 
not secure. In particular, the gas in the absorbers toward 
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Fig. 20. — log U vs. log 7V(HI) for the systems expected to be pre- 
dominantly photoionized. The dashed line indicates a theoretical 
prediction based on comparisons between CDM numerical simula- 
tions and the low redshift Lya forest (Dave & Tripp 2001). 



z = 0.36332 and z — 0.49510 may be in multiple ion- 
ization phases. In one case, we have only placed a lower 
limit to the U value; this is designated by the up arrow. 
Examining the figure, one notes a general decline in U 
with increasing AT (HI). Performing the conservative gen- 
eralized Spearman test, we calculate the significance to be 
85%. Similarly, a Pearson test adopting the lower limit as 
a value rules out the null hypothesis at 97% c.l. and the 
results are more significant if we exclude the 'anomalous' 
absorber at z = 0.36332. 

Overplotted in the figure (dashed-line) is a prediction for 
the correlation between U and AT (HI) at z = 0.2 derived in 
the following way. First, we adopt the Dave et al. (1999) 
relation between volume density nn and A(HI) for low z 
Lya absorbers: 

1 + 5 = — = —S- tx A(HI) 0,7 10~°' 4z (2) 
P n H 

where 

fi H = MVc(l + zf (3) 

is the mean hydrogen density of the universe. Second, 
we assume a constant HI ionizing flux. This gives U oc 
A(HI) -0 ' 7 and sets the slope of the dashed line in Fig- 
ure 20. Finally, the normalization of the line is determined 
by assuming a value for the HI photoionization rate; the 
central value reported by Dave & Tripp (2001) from their 
analysis of the statistical properties of the z < 1 Lya forest 
(Thi = 10 -13 - 3 photons s _1 for z = 0.17) gives the result 
shown here. 

The correspondence between observation and theory is 
impressive. The curve is not a fit to the observations, it is 
a a comparison of prediction from numerical simulations 
with observed properties of the low z Lya forest. It is es- 
pecially noteworthy that the curve results from an analysis 
of the low z Lya forest while our results are derived from 
the metal-line systems presented in this paper. Of course, 
the metal-line systems are simply a subset of the Lya for- 
est although the metal-line sample exhibits A(HI) values 
which exceed the principal range of the Dave et al. (1999) 
analysis. Furthermore, we caution that the majority of 
metal-line systems have AT (HI) « 10 14 ' 5 , i.e. the extremes 
of the parameter space described by Figure 20 are sparsely 
sampled. Nevertheless, the results presented in Figure 20 
support the model of the IGM inferred from cosmological 
simulations. 

Pushing the analysis one step further, we find a startling 
result regarding the total hydrogen column density A(H) 
of the low z Lya forest. Examining Figure 21 and column 
6 of Table 17, one notes that the predicted A(H) values for 
the metal-line systems based on photoionization modeling 
are nearly independent of AT (HI). For all but one of the 
absorbers 12 , logA(H) = 18.75 ± 0.3 dex. While this result 
could be biased by our reliance on metal-line systems (e.g. 
lower A(H) systems may not be detected unless the gas 
had a very low ionization fraction), the roughly constant 
A(H) value is a natural consequence of the correlation 
between uh and A (HI). Higher A (HI) implies higher nu 

12 The one exception is the unusual OVI absorber at z = 0.3633 
which shows a very weak Lya profile apparently offset from the 
metal-line transitions. We include it here, but note its peculiar na- 
ture. 
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which implies lower U which implies a lower ionization 
fraction and, finally, a lower A r (H)/7V(HI) ratio. In this 
manner, the correlation between nn and A^(HI) can lead 
to a roughly constant AT(H) value. 

The dashed line in Figure 21 presents the predicted 
N(H) values as a function of AT(HI) assuming (1) the 
EUVB model Q at z = 0.3 and (2) log U = -3-[JV(HT) /10 16 - 5 f 
with [3 — 0.7. The latter is equivalent to assuming a con- 
stant intensity for the EUVB radiation field and adopting 
the n H oc N(mf' 7 relation from Dave et al. (1999). We 
reach the unsettling conclusion that all Lya 'clouds ' have 
roughly the same AT(H) value. We stress that this result 
is sensitive to the value of /?; larger (0 > 1) or smaller 
((3 < 0.5) values imply steeper relations between 2V(H) and 
AT(HI). One seeks a physical interpretation - if it exists 
- for identifying N(H) » 10 18l5 ~ 19 cm" 2 as a fundamental 
scale for the low redshift Lya forest. This N(H) value can 
be identified with a comoving length scale by adopting the 
mean hydrogen density of the universe at z = 0. Taking 
n b = 0.04/if 2 (O'Meara et al. 2001), AT(H) = 10 19 cm- 2 
corresponds to i^ya ~ 17 Mpc at z = 0. The observed 
variations in A" (HI), therefore, result from the compres- 
sion of gas along this pathlength in overdense regions. 

It is worth noting that the description of the IGM de- 
veloped by Schaye (2001) based on Jeans length argu- 
ments also predicts a shallow dependence of AT(H) with 
nn and iV(HI). In a future paper, we will investigate this 
N(H) conspiracy through an analysis of numerical sim- 
ulations and analytic arguments. If the near constancy 
of iV(H) is confirmed observationally, it could challenge 
models of the Lya forest at AT(HI) < 10 16 cm -2 which 
associate the 'clouds' with individual galaxies (e.g. galac- 
tic halos; Lanzetta et al. 1995; Chen et al. 1998; Chen, 
Lanzetta, & Webb 2001; Manning 2003). 

In passing, we also note that because Dave et al. (1999) 
find fiH oc iV(HI) u ' to at least z = 3, one predicts the in- 
dividual Lya forest 'clouds' have a characteristic AT(H) at 




N(HI) 

Fig. 21. — Inferred total hydrogen column densities -/V(H) as a 
function of N (HI) value for the metal-line systems toward PKS0405— 
123 which we believe are predominantly photoionized. The dashed 
line in Figure 20 maps to the dashed line shown here. Remarkably, 
both observation and theory suggest the Lyo forest clouds have 
roughly uniform iV(H) value for iV(HI) = 10 12-16 cm -2 . 



any given redshift from z = — 3! The characteristic value 
of N(R) does vary with redshift owing to evolution in the 
EUVB intensity, the redshift dependence of Equation 2, 
and the increasing density of the universe (Equation 3). 
In general, AT(H) is predicted to increase from z = to 
z = 3. 

Before concluding this sub-section, we note that we searched 
the literature for additional photoionized, metal-line sys- 
tems to test the above conclusions. Unfortunately, this re- 
sulted in only two cases. One of these - the absorber at z = 
0.068 toward PG0953+415 (Savage et al. 2002, ; N(BI) 
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14.35 , 
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T.35) - has AT (HI) and U values which 



place it along the correlation expressed by the metal-line 
systems toward PKS0405-123 (Figure 20). Although the 
authors did not report an AT(H) value, adopting their ion- 
ization parameter we estimate AT(H) = 10 1 cm -2 , near 
cm~ 2 . The only other example is the absorber at 
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z = 0.0053 toward 3C 273 which lies within the Virgo su- 
percluster (Tripp et al. 2002). Taking their single-phase 
solution, this absorber falls below the relation in Figure 20, 
although its unique surroundings may imply an unusual 
physical environment (Stocke et al. 2004). It may also 
be noteworthy that the two outliers from the U, AT (HI) 
relation fall beneath the curve. This is consistent with the 
results of Dave et al. (1999) who emphasized the simu- 
lations show a significant scatter in the 1 + 6 oc AT (HI) 
correlation at low z with the majority of 'outliers' identi- 
fied above the relation (i.e. lower U). 

6.4. Chemical Enrichment of the Low z IGM 

An active area of current research is the study of the 
chemical enrichment of the IGM (e.g. Songaila 2001; Pet- 
tinietal. 2003; Schaye et al. 2003). These studies test sce- 
narios of Population III star formation (e.g. Gnedin & Os- 
triker 1997; Wasserburg & Qian 2000) as well as feedback 
processes in early galaxy formation (Aguirre et al. 2001). 
Preliminary studies with lOm-class telescopes argued for 
significant enrichment in Lya 'clouds' with Af (HI) > 10 14 cm" 
(Tytler et al. 1995) but a stacked spectrum of lower H I 
systems showed no corresponding C IV absorption (Lu et 
al. 1998). Recently, analyses of the 'pixel technique' have 
argued for metal enrichment to very low overdensity (e.g. 
Ellison et al. 2000). Recently, Schaye et al. (2003) have 
argued that the Lya forest has a median carbon metallicity 
[C/H] = -3.47t°;^ + 0.08^(2- 3) + 0.65T°;i°(log S- 
0.5) where S is the overdensity of the gas and is related 
to N(KI) as described in Equations 2,3. Simcoe, Sar- 
gent, & Rauch (2004) have reached similar conclusions, 
although with w 0.3 dex higher values for [O/H] and a 
less dramatic dependence on overdensity 13 . Very little 
research, in comparison, has been pursued at low red- 
shift. Buries & Tytler (1996) surveyed OVI absorbers 
at z ~ 1 and argue that the metallicity in these systems is 
[O/H] > -2.4 dex. Similarly, Barlow & Tytler (1998) ex- 
amined a large sample of C IV absorbers and found a mean 
metallicity [C/H] rs —1.5 dex. These results, imply higher 
enrichment levels within the Lya forest at low z and there- 
fore a significant evolution in the chemical enrichment of 
the IGM. 

13 These differences are primarily due to the fact that the authors' 
fiducial models assumed EUVB radiation fields with differing spec- 
tral shapes. 
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The analysis presented in this paper (summarized in Ta- 
ble 17) can be used to evaluate the chemical abundances 
of the low z Lya forest for comparison with high z. We 
will focus here on absorbers with logA(HI) > 14 which 
includes the absorbers without significant metal-line tran- 
sitions (§ 5). For the systems with detectable metal- line 
absorption, the median metallicity lies at « 1/10 solar en- 
richment if the absorbers are predominantly photoionized. 
A significant fraction of these absorbers exhibit chemical 
enrichment comparable to the Milky Way and its Magel- 
lanic Clouds. Another significant result is the large scatter 
in chemical enrichment among the absorbers. The systems 
span at least two decades in metallicity ranging from ap- 
proximately solar metallicity to < 1/100 solar. Although 
UV spectroscopy has poorer resolution and S/N than op- 
tical spectra, we emphasize that the relatively high metal- 
licities derived here arc not the result of observational 
bias. For most of the systems, these transitions are sat- 
urated or detected well above threshold. If this gas had 
10 x lower metallicity, several transitions would still be 
detectable, although constraints on the ionization state 
would be weaker. If we include the w 5 absorbers without 
significant metal-line absorption and assume their enrich- 
ment is below — 1.5dex, then the median metallicity lies 
closer to 1/30 solar. In either case, the median metallic- 
ity is consistent with the results presented by Barlow & 
Tytlcr (1998) based on their analysis of CIV absorption 
in the z ~ 0.5 Lya forest. 

Altogether, the results from the PKS0405-123 sightlinc 
indicate that Lya clouds with A(HI) > 10 14 cm~ 2 show a 
median value of 1/30 — 1/10 solar and a log normal dis- 
persion of w 1 dex. This median metallicity is significantly 
higher than the value ascribed to Lya 'clouds' with com- 
parable A(HI) values at high z (e.g Songaila 2001; Schaye 
et al. 2003). In fact, the value and scatter are remarkably 
similar to values derived for z w 2 damped Lya systems 
(e.g. Prochaska et al. 2003). Within the CDM paradigm, 
however, the 'clouds' in the low redshift Lya forest have 
much higher overdensity S compared to high z Lya clouds. 
If a correlation between metallicity and overdensity exists, 
we must be careful to compare mctallicities at the same 
overdensity. We can derive the metallicity of absorbers 
with comparable overdensity in the z ~ 3 universe by (i) 
calculating 1+5 for an absorber with AT (HI) = 10 145 cm~ 2 
using Equation 2 and (ii) evaluating the metallicity [C /H] 
at this overdensity using the results presented by Schaye 
et al. (2003). We find [C/H] = -2.9 at z = 3 with a log 
normal scatter of 0.55 dex. This metallicity and scatter 
are lower than the sample of absorbers analyzed in this 
paper. Even if we adopt non-solar relative abundances 
(e.g. [Si/C] = +0.7; Aguirre et al. 2004), the metallic- 
ity of these low z Lya absorbers greatly exceeds absorbers 
with comparable overdensity at high redshift. 

In short, our observations suggest a significant enrich- 
ment of the 5 w 20 IGM between z = 3 and the present 
universe. This is in stark contrast to current claims that 
the IGM has at most a moderate metal enrichment from 
z = 2 to 5 (Songaila 2001; Pettini et al. 2003; Schaye et 
al. 2003). Of course, our range of redshifts corresponds 
to w 10 Gyr wherein significant galactic enrichment has 
occurred. It apparently requires that the overdense IGM 
is polluted by processes related to galaxy formation, e.g., 



winds driven by supernovae and/or mergers (e.g. Cox et 
al. 2004). By inference, this weakens arguments that 
the IGM has been enriched primarily by an initial genera- 
tion of Population III stars (e.g. Gncdin & Ostriker 1997; 
Wasserburg & Qian 2000), although it docs not preclude 
the occurrence of both early and ongoing enrichment. 

Given a correlation exists between A (HI) and rijy, one 
might expect correlations between A (HI) and metallicity. 
For example, regions of higher overdensity will on average 
be closer to galaxies and therefore may be at higher metal- 
licity than less dense gas. As noted above, Schaye et al. 
(2003) report a correlation between metallicity and over- 
density when matching their simulations against C IV ab- 
sorption in the high z Lya forest. In the sample of z < 0.5 
metal-line systems toward PKS0405-123, however, there 
is no obvious correlation between A(HI) and metallicity. 
While the partial LLS at z = 0.16710 does show near solar 
abundance, the lower A(HI) absorbers at z = 0.49510 and 
0.36080 have comparable metallicity. Furthermore, there 
is no apparent correlation of metallicity with ionization pa- 
rameter or redshift. Because these conclusions are based 
on a small sample, an underlying trend may be masked 
by only a few outliers or systematic effects. A goal of fu- 
ture studies will be to obtain complete observations for a 
sample of A > 20 absorbers with a wide range of A(HI) 
values. 

In a smaller sample of the PKS0405-123 metal systems, 
we were able to constrain the N/O or N/Si relative abun- 
dance. Table 17 summarizes these results. The major- 
ity of absorbers are consistent with the general pattern of 
N/O values expressed by local measurements of HII re- 
gions and stars (Henry, Edmunds, & Koppen 2000). The 
obvious exception is the absorber at z = 0.49510 which 
shows a super-solar metallicity and an upper limit to N/O 
of —0.6 dex under the assumption of photoionization. The 
low N/O value is characteristic of the values measured in 
starbursting galaxies (e.g. Contini et al. 2002) and sug- 
gests this gas was enriched by an outflow from such a sys- 
tem. Future observations of N V and a deeper search for 
galaxies at z = 0.495 will help clarify this picture. 

7. SUMMARY 

We have presented an analysis of nine metal-line ab- 
sorption systems toward PKS0405-123 using the combina- 
tion of FUSE and HST/STIS ultraviolet spectroscopy. We 
derive H I column densities A(HI) from curve-of-growth 
analyses and examine ionic ratios to constrain the ioniza- 
tion mechanism of the gas. These ratios are compared 
against single-phase model predictions for collisional ion- 
ization and photoionization to constrain the temperature 
and/or ionization parameter U and thereby infer physical 
conditions of the absorber (density, metallicity ionization 
fraction) . 

We identify six O VI absorbers to a 3<7 equivalent width 
limit of 30mA along a non-contiguous redshift path Azo vi = 
0.38. This gives dN/dz\ovi — 16±6, consistent with other 
low z sightlincs. These systems exhibit metallicities [O/H] 
from —1.5 to Odex which roughly matches prediction for 
the warm- hot intergalactic medium (WHIM) . Half of the 
OVI absorbers, however, show transitions which suggest 
the gas is photoionized or has multiple ionization phases. 
This may argue against their membership in the WHIM, 
although a quantitative assessment awaits a careful as- 
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sessment of the ionization state of the gas in the numer- 
ical simulations. We also survey the WHIM through a 
search for NeVIII absorbers which should probe gas at 
temperatures T « 10 6 K. Over a non-contiguous path- 
length AznoVIH = 0.11 wc find no such systems and esti- 
mate dN/dz\^ e vni < 40 to an 3<t equivalent width limit 
of w 40mA. 

The photoionized metal-line systems show decreasing 
ionization parameter with increasing iV(HI) value. Both 
the slope and normalization of this correlation match a 
prediction derived from previous analysis of the Lya for- 
est with cosmological simulations (Dave & Tripp 2001). 
This correspondence lends support to paradigm of the 
IGM described by numerical simulations. Furthermore, 
the observed and predicted U vs. iV(HI) correlation im- 
ply a startling 'iV(H) conspiracy': the entire population of 
Lya clouds have roughly identical total hydrogen column 
density, i.e. independent of 7V(HI). We argue this conspir- 
acy holds at all redshifts although with the characteristic 
N(H) value increasing with redshift. 

The median metallicity of absorbers with AT(HI) > 10 14 cm 
is « 1/30 — 1/10 solar metallicity with a large (« 1 dex) 
dispersion. This metallicity is consistent with a previous 
analysis of the low z Lya forest and indicates a significant 
enrichment of the IGM since z ~ 2. There is no obvi- 
ous correlation between metallicity and any other physical 
property of the gas. Finally, we present a small sample 
of N/O measurements which are generally consistent with 
H II regions and stars of comparable metallicity. 
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